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Three-dimensional subsonic diffusers are relatively 
common in modern airbreathing propulsion systems. The shape 
of the duct cross-r action may vary in the streamwise 
direction, and offset bends ore often present. Strong 
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secondary flows can thus ho generated* In addition* since 
the flow is diffusing* the boundary layer may grow until its 
thicHness is comparable to the duct radius. These phenomena 
have important effects on overall diffuser performance and 
engino/inlet compatibility . 

Conventional boundary layer methods cannot properly 
analyse these comploK flows. A complete three-dimensional 
Navier-Stokes analysis could be used* but entremely large 
amounts of computer time and storage would be required for a 
detnilcil solution, thus precluding the use of such an 
analysis for routine calculation using present computer 
technology , 

However* if it is assumed that a primary flow direction 
con be identified, certain approximations can be made in the 
Havier-StoHos equations* resulting in a set of equations for 
fully viscous* subsonic* compressible flow that can be 
solved by forward marching. One such method, designated 
PEPSIG* was developed by Briley and McDonald^, and recently 
modified by Levy, McDonald, Briley, and Kreskovsky^. This 
paper presents the results of the first in a series of 
studies to evaluate the PEPSIG analysis. 

In this analysis, the flow is computed by a spatial 
marching procedure which solves an approximate form of the 
Navier-Stokes equations. Three basic assumptions are made. 
First, it is assumed that the flow is primarily in the 
direction of the duct centerline, with transverse secondary 
flow. Second, the pressure field is represented by the sum 
of a previously determined three-dimensional pressure field 
and a one-dimensional pressure correction computed as part 
of the marching analysis. And finally, second derivatives 
in the primary flow direction are assumed negligible. 
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The derivation of the governing equations actually 
solved in the PEPSIG analysis is fairly compleKi and has 
been presented elsewhere^. For completeness, however, a 
brief discussion, is presented here. The equations are first 
written in local Cartesian coordinates X2» and K3, with 
corresponding velocities Uj^, U2, and U3, The K3 direction 
is the primary flow direction, and is aligned with the duct 
centerline at each marching step. The xj^ and X2 directions 
define the transverse, or secondary flow plane. The 
equations are then transformed into general non-orthogonal, 
body-fitted, curvilinear coordinates y^^, y2, and y3- 

The X3 momentum equation is then given, using tensor 
notation, by: 
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Here P3_p “ P3_p (xj^.x^jx^) is the known three-dimensional 
pressure field and “ P^-d ^^3) correction compute 

during the marching solution. In the present study, p3„p was 
computed using a three-dimensional potential flow analysis . 


The differential continuity equation is written in 
terms of a scalar potential <|> as: 
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The and K 2 momentum equations are combined to form 
on equation for vorticity transport in the K 3 direction and 
a cross -plane stream function equation. After some 
simplifications, the vorticity equation is given by: 

“J ■ P 3^2 < 

and the stream function equation by: 
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The energy equation is eliminated by assuming constant 
total enthalpy. A mixing length turbulence model is used, 
uith the distribution of mixing length given by the Buleev 
formula^. The boundary conditions used result in no-slip at 
the walls, and symmetry at the symmetry plane. 

The above equations are solved by forward marching from 
an initial station, where the flow profiles are known, using 
an implicit numerical technique. Details of the solution 
procedure are given in references 1 and 2 . 


COMPUTED RESULTS 


Circular S-Duct 

As an initial test of the ability of the PEPSIG 
analysis to correctly predict the basic physical phenomena 
present in curved ducts, a simple circular S-duct was run. 
The duct had a diameter of .61 m. (2 ft.) and the centerline 
was offset one diameter in a streamwise distance of five 
diameters. The flow was laminar, with an entrance Mach 
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number of »2 and a Reynolds number based on diameter of 
200t). The initial boundary layer thich.ness uas 3.0 cm. (.1 
ft,,). This case was run uith a 20 h 20 mesh in the 
cross-plane, and streamuise stations. The calculation 
required 16 minutes on on IBM 570/3033. 

Computed total pressure contours are shown at several 

stations along the duct in figure 1. The boundary layer 

growth along the top of the duct in the first bend can bo 

clearly seen, In the second half of the duct, this region 

of low energy flow espands greatly, forming a pocket of low 

total pressure at the top of the duct. This phenomenon has 

4-5 

been seen in several enperiments . A comparison between the 
analytical results and the results from one of these 

A 

experiments is presented in figure 2. Here total pressure 
contours at the e.uit plane are plotted for both coses. 

Although only a qualitative comparison can be made since the 
tested case was for a different geometry and flow conditions 
than the computed case, the basic physics of the flow appear 
to have been computed correctly. 

The mechanism by which the pocket of low energy flow is 

A 

formed has been presented by Bansod and Bradshaw . This 
mechanism is present in the analytical results, and can be 
explained by examining the cross-flow velocity vectors, as 
shown in figure 3. Except for the results at station 5, 

Which is 1/A of the way through the first bend, the vector 
plots are shown at the same stations as the total pressure 
plots in figure 1. In the first half of the duct (stations 
1-21), the classic secondary flow pattern for flow through a 
curved pipe is set up. The low speed boundary layer fluid 
migrates circumferentially from the pressure surface toward 
the suction surface, and the essentially inviscid core flow 
responds to centrifugal forces. A pair of contra-rotating 
Vortices is thus established. In the second bend (stations 
21-Al), these vortices are intensified, and become centered 
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near the toi> of the duct. As the lou energy fluid 
npproachos the symmetry plane from both sides r it erupts 
from the surface forming the pocket of lou total pressure 
characteristic of S-duct flows, 

To evaluate the ability of the analysis to compute the 
flow thrc:'‘!h a realistic diffuser configuration! 
calculations upre made for subsonic flow in the F'“16 inlet'^, 
A schematic diagram of this inlet is given in figure 
Cross-sections are shown at various stations along the 
inlet. The station numbers are fuselage stations from 
reference 6. The inlet has a generally S-shaped centerline* 
with the first and second bonds separated by a straight 
section. The total offset is .616 m, (2.02 ft.) and the 
length from the throat to the engine face is ^.7 A m. (15,55 
ft,). The cross-section shape transitions from slightly 
smile-shaped at the throat to circular at the engine face. 
The overall area ratio is 1.33. 

The non-orthogonal, body fitted mesh generator 
currently used in the PEPSIG analysis requires 
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super-elliptic cross-section shapes . Meshes can thus be 
generated for ducts witli cross-sections that vary from 
nearly rectangular or square (with a large super-ellipse 
exponent) to elliptical or circular (with an exponent of 2). 
The actual F-16 cross-sections could therefore not be 
modeled exactly. Instead, the following procedure was used. 
At each streamwise station, the major axis of the 
super-ellipse was chosen to match the maximum spanwise 
dimension of the actual cross-section. The minor axis used 
was the gapwise distance in the symmetry plane of the inlet. 
The super-ellipse exponent was then computed to give the 
correct area, but was not allowed to fall below 2, 
corresponding to an ellipse. The result was the baseline 
configuration shown in figure 5a. 
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A calcuXation was made for turbulent flow in this 
configuration at a Reynolds number per meter of I.AAhIO^ 

per foot) and a free stream Mach number of .9, The 
initial boundary layer thickness specified was .91 cm. (.36 
in,). These flow conditions correspond to one of the test 
conditions of reference 6. A 2 <)k 2A mesh was used in the 
cross-plane I with 250 streamwise stations. The calculation 
required about 2 1/2 hours on an IBM 370/3033. 

Computed total pressure contours are shown at several 
stations along the duct in figure 6. Boundary layer growth 
around the entire circumference is evident through the first 
half of the duct, About halfway through the duct, the 
boundary layers at the top and bottom start growing more 
rapidly, while near the sides the boundary layer becomes 
thinner . 

The surprising aspect of this calculation was the 
degree of symmetry between the contours in the top and 
bottom halves of the duct. Based on the S-duct calculations 
presented earlier, as well as enperimental results for 
S-ducts, it had been expected that secondary flow vortices 
would be generated by the centerline curvature and would 
cause a shift in the total pressure contours. A comparison 
of computed and experimental exit plane total pressure 
contours is presented in figure 7 , The data are from 
reference 8, and ore for a somewhat different geometry and 
incoming flow than was used in the analysis. (In the test, 
the inlet had an external ramp forebody and a compressor 
face hub, and was run at a 1 degree angle of attack. The 
Reynolds number was A. 9x10 per meter (1.5x10° per foot) and 
the free stream Mach number was .8). The experimental 
results, however, are similar to the computed results in the 
degree of symmetry between the upper and lower halves of the 
inlet. It should be noted that the computed boundary layer 
grows very quickly during the first few marching steps. 
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This is believed to bo caused by appreciable numerical 
diffusion in this thin boundary layer calculation as a 
result of inadequate mesh “^solution. This resulted in 
larger calculated total pressure losses than were measured 
enporimentally. With proper mesh resolution » tias problem 
should be alleviated* 

It was then hypothesised that the effect of the change 
in cross-section shape was much more important in this flow 
than the effect of the centerline curvature* To verify this 
hypothesis, two additional geometries were run. The first, 
shown in figure S5b, had the same distribution of 
cross-section shape as the baseline configvirntion, but with 
a straight centerline* The second, shown in figure 5c, had 
the same centerline shape and area distribution as the 
baseline configuration, but with circular cross-sections. 
Those configurations were chosen to isolate the effects of 
shape change and centsrli.ns curvature. The flow conditions 
and grid were the same as for the baseline calculation. 

The total pressure contours for the straight centerline 
configuration are shown in figure 8, and are virtually 
identical to those for the baseline configuration. The 
results for the circular cross-section case, shown in figure 
9, do show a slight shift in the total pressure contours, 
indicating that the centerline curvature is strong enough to 
generate at least a weak secondary flow vortex* This is 
confirmed when the cross-flow velocities are examined. The 
computed secondary flow velocity vectors at the exit plane 
for all three F-16 configurations are shown in figure 10. 

The secondary flow patterns for the baseline and straight 
centerline configurations are basically similar. Although 
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some influence of centerline curvature can be seen in the 
baseline results, no vortex pattern is set up. The circular 
cross-section configuration, however, has an obvious 
secondary flow vortex at the exit plane. 
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Additional insight can he gained by eKamining the 
circumferential pressure gradient in the baseline F-16 
configuration. In figure 11a the inviscid static pressure 
is plotted for a cross-section at about the mid-point of the 
duct. The pressure is higher in the region of the 
converging sides of the cross-section than at the diverging 
top and bottom. The low energy boundary layer fluid is thus 
driven in hO-th directions away from the sides of the inlet, 
causing the boundary layer thickness to decrease there and 
to increase at the top and bottom of the inlet. By the time 
the second bond is reached, most of the cross-section 
transitioning has been completed. The circumferential 
pressure distribution in the second bend, shown in figure 
11b, is similar to that found in the second bend of a 
constant area S-duct. However, the boundary layer has 
become so thin near the converging sides of the duct that 
the boundary layers in the top and bottom halves of the 
inlet are essentially isolated from each other, low energy 
boundary layer fluid cannot migrate through this thin 
boundary layer region, and therefore no secondary flow 
vorteK is formed. In the circular cross-section 
configuration, of course, the boundary layer does not thin 
out at the sides of the duct, aud a vorteK pattern is 
generated. Thus, although the centerline curvature in the 
baseline F-16 configuration is strong enough to produce a 
secondary flow vortex, the changing cross-section shape 
suppresses its formation. 

GONCLUniNG PEHARkS 

1 , The PEPSIG analysis gives at least qualitatively 
correct predictions of the development of secondary flow 
vortices and total pressure distortion in S-shaped ducts. 

2, More detailed quantitative comparisons with 
experiment are needed to determine the accuracy of the 
analysis . 
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5. Mesh studies ore needed to determine the resolution 
required for thin boundary layer calculations. 

4. The analysis can fee u»ad to help understand the 
development of secondary flotj phenomena in compleK 
configurations. 

5. There is a critical need for detailed enperimental 
data on secondary flow d^elopmont in ducts with centerline 
curvature and/or cross-section transitioning for proper code 
evaluation. 
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Figure 9. - Computed totel pressure contours (or F-16 Inlet, circular cross-section 
configuration, 



SELINE STRAIGHT CENTERLINE CIRCULAR CROSS-SECTION 

Figure 10. - Computed secondary flow at exit plane o( F-16 Inlet configurations. 





